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A  STUDY  OF  THE  VIBRATIONS  OF  SHALLOW  SPHERICAL  SHELLS 

ABSTRACT 

In  this  report  there  is  presented  a  comparison  of  the  theoretically  and 
experimentally  determined  frequencies  and  zero  displacement  lines  for  the 
normal  mode  vibrations  of  shallow  spherical  shells.  Results  are  given  for 
the  clamped  edge  and  the  momentless  edge  boundary  conditions. 

Calculations  of  frequencies  and  loci  of  zero  displacement  components  for 
various  shell  sizes  were  made  on  an  IBM  650  computer  for  the  rotationally 
symmetric  vibrations  only,  as  the  theory  is  not  available  for  the  case  of 
asymmetric  vibrations. 

Experiments  were  performed  to  determine  the  vibration  characteristics 
both  for  the  symmetric  and  the  asymmetric  vibrations.  Several  model  shells 
constructed  of  steel  were  used  for  the  purpose. 

Results  of  the  investigation  definitely  demonstrate  the  value  of  the 
theory  for  the  determination  of  the  dynamic  response  of  shallow  shells. 
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INTRODUCTION 


1* 

In  a  previous  paper  ,  the  equations  of  motion  for  flexural  vibrations 

of  orthotropic  shells  were  developed  and  simplified  so  as  to  obtain  solutions 

for  the  case  of  shallow  shells.  In  particular,  it  was  shown  in  that  paper  how 

the  equations  reduce  to  cover  the  case  of  isotropic  shells  for  which  some 

computed  frequencies  were  given  for  both  the  shell  with  clamped  edge  and  with 

momentless  edge.  The  equivalent  equations  for  the  clamped  isotropic  shell 

2 

were  given  previously  in  a  paper  by  Eric  Reissner 

Because  of  the  importance  of  the  problem  and  because  some  doubt  exists 
at  the  present  time  as  to  the  applicability  of  the  theoretical  results,  ex¬ 
periments  on  shell  models  were  considered  to  be  needed.  In  addition  a  com- 

3 

putational  program  was  also  needed  to  check  the  previous  calculations  and 
to  extend  them  to  cover  a  larger  number  of  shell  parameters. 

The  present  report  presents  experimentally  determined  frequencies  and 
zero  displacement  loci  for  both  symmetric  and  asymmetric  vibrations.  It  also 
provides  theoretically  determined  frequencies  and  zero  displacement  loci 
computed  with  an  IBM  650  computer  for  the  symmetric  vibrations. 

THEORY 

The  thin  shell  theory  for  small  amplitude  vibrations  which  is  under 
consideration  in  this  report  has  been  amply  discussed  in  texts  by  A.  E.  H.  Love, 
S.  Timoshenko,  and  others.  It  is  basically  a  stress  resultant  theory  in  which 
the  resultant  of  the  stress  distribution  through  the  thickness  of  the  shell 
at  any  point  may  be  represented  as  a  force  and  a  moment  per  unit  length  of  the 
line  formed  by  the  intersection  of  the  middle  surface  and  a  plane  normal,  to 
that  surface.  In  such  a  case  the  defining  differential  equations  of  motion 
are  linear.  Despite  this  fact  however  they  are  in  general  quite  difficult  to 
solve  and  consequently  very  few  solutions  to  the  dynamical  equations  for  shells 
have  been  published.  It  turns  out  that  in  the  important  case  of  shallow  shells, 
the  equations  can  be  solved  and  the  frequencies,  as  well  as  lines  of  zero 


Superscript  numbers  refer  to  references  at  end  of  report. 
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displacement,  can  be  found.  The  frequency  and  displacement,  equations  have  been 
given  in  a  previous  paper1.  A  summary  of  equations  are  provided  in  the  Appendix 
to  the  present  paper. 

In  the  previous  paper,  frequencies  computed  with  the  aid  of  electric  desk 

4 

calculators  and  Tables  of  Functions  were  tabulated.  There  the  assumed  shells 
had  chord  length  of  12,  thickness  of  l/4,  and  sagittas  of  0.5,  1.0,  and  1.6. 

In  the  present  study  frequencies  and  zero  displacement  lines  of  the  normal 
modes  were  obtained  with  the  IBM  650  computer  and  Tables  of  Functions  ^or 

shells  of  chord  length  12,  thicknesses  of  l/l6,  l/8,  l/4  anrt  sagittas  of  0.5 
and  1.6.  The  calculations  were  carried  out  for  clamped  and  momentless  edge 
conditions.  All  of  these  calculated  results  are  given  in  the  Appendix  in  Tables 
II-IV. 

VIBRATION  EXPERIMENTS  WITH  THIN  SHALLOW  SPHERICAL  SHELLS 

For  the  purpose  of  making  experiments  two  model  shells  were  designed  and 
constructed.  An  arrangement  of  the  shell  with  the  device  used  to  provide  a 
clamped  edge  boundary  is  shown  in  Fig.  1.  The  momentless  edge  boundary  condi¬ 
tion  was  obtained  by  using  the  clamping  ring  the  same  as  for  the  clamped  case 
but  machining  a  deep  peripheral  groove  0.84  times  the  thickness  of  the  shell 

and  very  close  to  the  clamping  ring,  in  the  same  manner  described  elsewhere  for 
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use  with  vibrating  plates  . 

The  shells  were  of  steel  and  they  were  vibrated  by  means  of  an  electro¬ 
magnetic  oscillator  shown  schematically  in  Fig.  1.  The  resonant  tuned  magnetic 

circuit  and  power  pack  to  obtain  sufficient  driving  power  had  been  previously 

0 

developed  .  The  lines  of  zero  displacement  were  found  with  the  use  of  a  crys¬ 
tal  pickup  whose  output  was  fed  to  two  of  the  plates  of  a  cathode  ray  oscil¬ 
loscope.  The  driving  frequeucy  from  an  audio -oscillator  were  fed  into  the 
other  two  plates  of  the  C.R.O.  The  technique  is  quite  satisfactory.  Its  use 
with  a  capacitive  type  pickup  instead  of  the  crystal  is  fully  described  in 
another  paper 

THEORETICAL  AND  EXPERIMENTAL  RESULTS 

The  experimentally  determined  frequencies  for  the  case  of  symmetric  vibra¬ 
tions,  both  for  the  clamped  edge  and  momentless  edge  boundary  conditions,  to¬ 
gether  with  the  corresponding  theoretical  results  are  given  in  Table  I.  In  that 
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Fig 


MOMENTLESS  EDGE _  CLAMPED  EDGE 


TABLE  I 


SYMMETRICAL  MODES  OP  FLEXURAL  VIBRATIONS 
OP  SHALLOW  ISOTROPIC  SPHERICAL  SHELLS 


FREQUENCIES  x(lO~3  rad. /sec. ) 


h 

s  ■  0 

S  = 

0.5 

S  =  1.0 

S  - 

1.6 

(in. ) 

n 

Theor. 

Theor. 

Exper. 

E9m 

Bm 

Exper. 

0.0625 

H 

1.1 

6.6 

11.8 

Q 

4.3 

8.8 

14.2 

20.0 

a 

9.6 

11.5 

17.8 

24.0 

0 

17.0 

18.0 

21.3 

0.125 

H 

2.2 

7.9 

6.5 

13.2 

19.1 

15.4 

a 

6.6 

11.0 

10.0 

17.6 

24.6 

22.3 

| 

19.1 

20.0 

17.1 

23.0 

29.9 

29.2 

B 

33.6 

34.4 

29.3 

36.0 

39.0 

42.2 

0.250 

B 

4.4 

9.1 

7.2 

15.7 

22.4 

18.5 

17.1 

16.4 

15.7 

22.0 

29.4 

22.3 

38.2 

36.6 

33.2 

40.0 

42.6 

34.2 

1 

66.0 

68.1 

56.5 

69.O 

70.4 

55.0 

0.0625 

1 

0.5 

6.4 

11.5 

2 

3.2 

6.8 

13.6 

19.4 

3 

7.9 

10.3 

17.7 

22.4 

4 

14.9 

15.9 

19.8 

0.125 

B 

1.1 

7.7 

6.2 

12.6 

18.7 

15.2 

B 

6.4 

9.8 

10.1 

17.6 

23.4 

19.6 

B 

15.9 

17.0 

16.5 

20.5 

29.7 

29.2 

B 

29.7 

26.7 

28.0 

31.8 

35.4 

40.6 

0.250 

B 

2.1 

8.8 

7.2 

15.5 

21.0 

17.3 

12.8 

14.3 

14.5 

19.6 

29.2 

21.0 

m 

31.8- 

32.4 

30.2 

34.0 

37.4 

33.7 

D 

59.2 

59.6 

51.6 

60.6 

61.6 

52.5 
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Table  h  is  the  thickness,  S  is  the  sagitta  length,  and  n  is  the  mode  number. 

The  experimentally  determined  loci  of  zero  radial  displacements  for  two 
thicknesses  of  shell  are  given  in  Tables  II  and  III  along  with  corresponding 
theoretical  results.  Definitions  of  symbols  are  in  Fig.  2. 

The  experimentally  determined  frequencies  and  loci  of  zero  displacement 
components  for  asymmetric  vibrations  are  given  in  Table  IV.  Because  the  loci 
for  these  vibrations  are  not  simply  circles,  sketches  of  the  loci  for  the 
several  modes  which  were  investigated  are  shown. 

DISCUSSION  AND  CONCLUSIONS 

It  is  concluded  that  the  results  of  the  investigation  demonstrate  that 
the  theory  describes  very  well  the  nature  of  the  vibrations  of  shallow  spher¬ 
ical  shells.  It  is  also  considered  that  the  same  theoretical  results  can  be 
used  for  describing  the  vibrations  of  any  very  shallow  shells  of  revolution 
which  have  positive  curvature  everywhere  and  which,  further,  have  the  same 
base  diameter  and  sagitta.  In  these  cases  the  shell  is  then  really  assumed 
to  be  equivalent  to  the  osculating  sphere  at  the  apex.  The  modal  shapes  for 
the  shallow  shells  are  similar  to  those  for  flat  circular  plates  but  the  fre¬ 
quencies  are  much  higher,  indicating  greater  stiffness. 

It  is  interesting  to  note  that  for  the  fundamental  vibration  of  the 
symmetric  type  there  are  two  nodal  circles,  one  at  the  boundary  and  another 
closer  to  the  center.  As  the  sagitta  or  rise  of  the  shell  approaches  zero 
the  latter  circular  node  approaches  closer  to  the  boundary  and  in  the  limit, 
which  is  a  flat  plate,  it  becomes  coincidental  with  the  boundary.  It  is  rather 
surprising  that  the  frequencies  for  the  clamped  edge  condition  and  the  moment¬ 
less  edge  condition  are  so  nearly  equal  in  the  various  modes.  Further  study 
should  be  made  of  this  point  in  order  to  assess  the  influence  of  membrane 
stresses  as  opposed  to  that  of  flexural  stresses.  It  would  be  pertinent  in 
such  a  study  to  compare  the  potential  energy  corresponding  to  stretch  and  the 
potential  energy  corresponding  to  flexure  as  was  done  by  Arnold  and  Warburton 

for  cylindrical  shells  in  order  to  assay  the  anomaly  of  complexity  of  nodal 
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pattern  as  related  to  frequency  of  vibration  . 
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TABLE  II 


LOCI  OP  RADIAL  ZERO  DISPLACEMENT  (SYMMETRIC  CASE) 


MOMENTLESS  EDGE  (h  =  0.125  in.) 


S  n  fl  theorv^-^^  T|  exper. 

(in.)  _ _n _ ^"-01 _ 


TABLE  III 


LOCI  OF  RADIAL  ZERO  DISPLACEMENT  (SYMMETRIC  CASE) 


~  _  ^i  Location  of  radial 
1  ”  zero  displacement  circle 
r0 

i  =  1,2,3  ••• 

wn  =  radial  displacement 

vn  =  tangential  displacement 

Fig.  2 
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TABLE  IV 
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APPENDIX 


COMPUTATIONS  FOR  DETERMINATION  OF  FREQUENCIES 
AND  NODAL  CIRCLES  OF  SYMMETRIC  VIBRATIONS 


To  carry  out  computations  the  following  dimensionless  parame¬ 
ters  were  introduced 

*n  ‘  »n  rof^S 
€  =  2f3  rQ/h 

(1) 

6  =  4fT  s/h 

n  =  x*/r0 


Applying  these  definitions  the  resulting  eigenvalue  problem 
is  defined  by  the  following  equations: 


1 

H 

1 

-v2K2  j  {x4-  (1- 

v2)i2  e2  +  262(1+V 

)}  -  (1+v)2  B2  x2  -  0  (2) 

CASE  A. 

Clamped  Edge. 

dw„ 

, ,  n 
w  =  -  » 

n  d”H 

v  =  0  at 

n 

ii 

«-• 

(3) 

Vxi) 

X1 

Jq(x2^ 

*2 

Vx3} 

x3 

Do  * 

Jx(x2) 

Jl(x3} 

o 

i 

J1(x1) 

Jll^) 

Jx(x3) 

42(l-vg)-x| 

ix(l-v2)-Xj 
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CASE  B.  Momentless  Edge. 

M .  -  0  at  T|  ■  1 


W  mV 

wn  n 


(5) 


M 


J0^X1^ 


JQ^X3^ 


J1(x1)(l-v)-x1J0(x1)  ,  J1(*a)(l-v)-*2,0(xa)  '  JlU3)(1'V)'X3J0(x3) 


Ji  [*iL_ 
e2(i-v2)-^ 


J^Xg) 

£2(l-v2)-x| 


Jx(x3) 

S2(l-v2)-x| 


-  0 


(6) 


The  eigen  functions  are  given  by 

wn  -  -  K  ^  +  I  W  *  J0<"*3>} 


(7) 


vn  = 


r  tMnx.) 

-(1+V)K  + 


J.  (TtXg)  Jn  (Tlxq) 

— im - £ - *  + 


?  +  U-v2)^  T  ’3 

(8) 


} 


where 

CASE  A.  Clamped  Edge, 


B 

B 


•2  A  -  A  . 

4  -  xf  u-v2)4t 


,2^2 


J^(Xg) 


(9) 


B 

B 


x2  -  x2 .  d-v 2 ln2  -  4}  Ji(xi} 

4  -  X2  ‘  (!-»■)«“  - 


(10) 
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CASE  B.  Moment less  Edge 


<l-v2)S2  -  x2 

X1J0^X2) 

h. 

(1-v2)^  -  x2 

Jl<xl> 

x2J0(Xi) 

(l-v2)*2  -  x2 

Jx(x3) 

x1J0(x3) 

(1-v2)52  -  X2 

x3^0^xl^ 

B1  _  /  xlJ0^x2^ 

+ 

B2  ‘‘X2Jq(x1) 

x3J0(xi)  J 

(ID 


All  computations  were  carried  out  on  an  IBM  650  computer  using  floating 
point  arithmetic.  The  necessary  exponential  subroutine  was  included  in  the 
floating  point  subroutine.  The  Bessel  function  subroutine  used  checked  to 
within  six  significant  figures  through  x  =  13  with  the  tables  listed  in  the 
references. 

The  values  of  parameters  used  were  determined  to  allow  direct  comparison 
with  experimental  results  using  shells  of  thickness  O.25  and  0.125  inches  with 
sagittos  of  0.5/  1.0  end  1.6  inches.  Table  I  in  the  Appendix  shows  the  cor¬ 
responding  values  of  e  and  8.  In  addition,  for  comparison  with  experiment 

2 

p  r  0  „  _io  2 

Tq  =  6  inches  and  — g —  *  8.819  x  10  sec  . 


The  computational  procedure  used  was  as  follows: 

1. 


2. 

3. 


Choose  e  and  8. 

^  2 
ChooBe  5  . 


Compute  the  six  roots  of  Eq.  (2)  calling  those  with  positive 
sign  x^,  x2  and  Xy  For  the  purposes  of  this  work  x^  was 
always  real  and  was  always  pure  Imaginary. 


4.  Using  the  chosen 


D  and  D„. 

C  M 


and  x^,  x2  and  x^  compute  the  values  of 


5.  Choose  another  value  of  $  and  repeat  steps  3  and  4. 


TABLE  I 


h 

e 

0.0625 

0.125 

0.250 

332.55 

I06.27 

83.13 

h  =  0.0625 

h  0  0.125 

h  =  0.25 

S 

5 

8 

5 

m 

55.42 

27.71 

13.85 

110.85 

55.42 

27.71 

177.36 

88.68 

44.34 
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The  zerou  of  the  resulting  graphs  of  and  vs.  i  yield  the  eigenvalues 

I  .  The  values  of  |  are  given  In  Tables  II  through  IV.  The  results  obtained 

&  u 

are  accurate  to  within  two  In  the  third  digit. 


In  the  eigen  function  confutation  the  functions 


w’ 

n 


w 

J&  ,  v* 

V  n 


(l+VjBj. 


were  computed  for  the  clamped  case  and 


V*  a  -  —  .  V*  =  - 

n  *  n 


(l+v)B„ 


(12) 


(13) 


were  computed  for  the  nomentless  case.  Zeros  of  the  eigen  functions  are  given 
in  Appendix,  Tables  II -IV. 


LOCI  OP  ZERO  DISPLACEMENT  COMPONENTS  (SYMMETRIC  CASE) 


|. 


2^ 


25 


ptwifc. 
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